Viable entanglement detection of unknown mixed states in low dimensions 



Thiago 0. Maciel 1 '* and Reinaldo 0. Vianna 1, t 

1 Departamento de Fisica - ICEx - Universidade Federal de Minas Gerais, 
Av. Antonio Carlos 6627 - Belo Homzonte - MG - Brazil - 31270-901. 
(Dated: July 6, 2009) 

We explore procedures to detect entanglement of unknown mixed states, which can be experimen- 
tally viable. The heart of the method is a hierarchy of simple feasibility problems, which provides 
sufficient conditions to entanglement. Our numerical investigations indicate that the entanglement 
is detected with a cost which is much lower than full state tomography. The procedure is applicable 
to both free and bound entanglement, and involves only single copy measurements. 
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I. INTRODUCTION 

Entanglement is in the heart of Quantum Informa- 
tion, and allows for promising spectacular applications 
like Quantum Computation and Unconditionally Secure 
Quantum Cryptography. Our understanding of entan- 
glement has largely grown in the last few years, but the 
experimental detection is still a daunting challenge. A 
recent and thorough review of the subject is provided by 
the Horodecki family [1]. 

Theoretically, the tool of choice to detect entanglement 
is an Entanglement Witness (EW) [2]. It consists of a 
Hermitian operator (W) with non-negative expectation 
values for all the separable states, but which can have 
a negative expectation value for an entangled state, in 
this case, the state is said to be detected by the EW. 
When it comes to experiments, EWs are not that good, 
for each state has its own optimal witness. The con- 
struction of the optimal EW depends on the knowledge 
of the state (see [3], for example). There does not exist 
an EW which detects all the entangled density operators 
acting on a given Hilbert space. Nevertheless, an EW 
can detect many states on a certain region of the state 
space, though it will be optimal just for a restricted fam- 
ily of states. Therefore, when some information about 
the state is known, an EW can be implemented. For ex- 
amples of experimental implementations of EWs, consult 
the references in [1]. 

Exploring collective measurements to estimate nonlin- 
ear functionals of quantum states, Walborn et al. [4] 
have experimentally measured the concurrence of un- 
known pure two qubit states, using two copies of the 
objective state. It has also been extended to the estima- 
tion of the concurrence of mixed states, and implemented 
experimentally [5]. In the case of rank-2 two-qubit states, 
it is also possible to measure the concurrence exactly by 
means of collective measurements on four copies [6]. 

Here we investigate measurements on single copies of 
unknown mixed quantum states. The method we propose 
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is simple, and shows to be effective in low dimensions. 
We are advocating that using sophisticated mathemat- 
ical tools to characterize entanglement in a data post- 
processing fashion, while keeping the experiment as sim- 
ple as possible, is efficient. Therefore it is an approach 
in the opposite direction of works like those in [4, 5], 
for example, where entanglement is directly measured in 
a very elaborate experiment. We have performed nu- 
merical tests in systems of two qubits, one qubit and 
one qutrit and two qutrits. In the case of two qutrits, 
we have investigated both bound and free entanglement. 
Though we discuss just bipartite entanglement, the for- 
malism can be straightforwardly applied to the multipar- 
tite case. The basic idea is to consider the state written 
in an orthonormal basis, which can be thought of as a 
generalized Block representation. Then the projections, 
which are the components of the generalized Block vector, 
are gradually measured. For each set of measurements, 
it is checked if there is enough information to infer en- 
tanglement. In the case of states with Negative Partial 
Transpose (NPT) , we check if there is no state with Pos- 
itive Partial Transpose (PPT) compatible with the mea- 
surements. In the general case, including entangled PPT 
states [7], we build an entanglement witness compatible 
with the measurements, using the method described in 
[3] . In the spirit of data post-processing, we mention the 
recent techniques independently introduced by Eisert et 
al. [8] and Guhne et al. [9] , that yield bounds to certain 
entanglement quantifiers, based on the measurement of 
non-optimal EWs. We note also that Badziag et al. [10] 
and Hassan et al. [11] have introduced interesting en- 
tanglement criteria based on the norm of the generalized 
Bloch vector. In the context of Quantum Key Distribu- 
tion, Curty et al. [12] have introduced a method to check 
the presence of entanglement by means of EWs built with 
previous measured data. 

In the next section, we introduce and illustrate our 
method. Section III discusses a possible choice of infor- 
mationanlly complete set of observables, and illustrates 
how our method would perform with projective measure- 
ments. In particular, that section makes clear that, com- 
pared to two qubits and two qutrits, it is not obvious how 
to choose the minimal amount of (or optimal) measure- 
ments in a 2 ® 3 system. Section IV offers more ques- 
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tions than answers. There we analyze the limitations of 
the method, studying three representative states, namely, 
one highly entangled and tow very low entangled states. 
As presented in this work, our method offers a yes/no an- 
swer about the entanglement of a unknown mixed state, 
but the calculations of section IV suggest that a further 
development of it could yield a good quantitative estima- 
tor of entanglement. Section V concludes the paper. 



II. THEORY 

Given a state represented by the density operator p, we 
want to check if it is entangled, without performing a full 
tomography. As matter of fact, we want to make the least 
possible number of assumptions about the state. We will 
present a strategy based on acquisition of partial informa- 
tion about the state, followed by data (post-)processing 
in form of Semi-Definite Programs (SDP). SDPs can be 
efficiently solved [13], and have exact solutions. As the 
numerical tests will show, it is effective in low dimen- 
sions. We focus on bipartite states in order to simplify 
the discussion, but the generalization of the formalism to 
multipartite states is straightforward. 

A state p, acting on the Hilbert space Hd = C da x C db , 
where d = d a x dj,, can be written as: 
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where the Pij are observables forming a complete basis 
in the Hilbert-Schmidt space, and = Tr(pPij) € 1Z. 
One possible choice for these observables is = 
being SU(d s ) matrices, i.e., gen- 
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eralizations of the Pauli matrices, where a ds stands for 
the identity matrix, and tqq = l/d. In this case, the 
state can also be written with explicit local and non-local 
parts, and we have an expression that can be thought of 
as a generalized Block representation, namely, 

P = ^(Id a ® Id b + r a ■ o d * <8> I db + Id„ <8> n ■ a d " + 

d\-\d\-\ 

»=i j=i 

where Id s is the d s x d s identity matrix, a ds are the ma- 
trices for SU(d s ), r s € IZ 4 ^ 1 , and finally tij G 1Z. Note 
that f a and r*b are the local parameters, defining the re- 
duced density matrices, namely: 



p a = Tr h p=— (I a + r-<7 d -) 



(3) 



where Tr^ is the partial trace on subsystem b, and an 
analogous expression for pi,. The non-local parameters, 



form a real matrix T, and are responsible for the classical 
and quantum correlations in p. Note that the parameters 
in Eq. 1 or Eq. 2 must be real in order to p to be 
Hermitian, but it does not guarantee its positivity. 

We will introduce some procedures to check the entan- 
glement of an unknown state, based on partial informa- 
tion about it. This partial information consists of the 
knowledge of some of the r^ (viz. Eq. 1), eventually 
enriched with some further characteristic of the state, as 
the fact that it is NPT, or its marginals are known (viz. 
Eq. 3). _ 

As it is well known, it is harder to check the entangle- 
ment of a bound entangled state than a free entangled 
one. From the theoretical point of view, it is easy to 
know if the latter are entangled, for they have negative 
partial transpose, which is known as the Peres-Horodecki 
criterion [2, 14]. But if the state is PPT, we need an 
entanglement witness. The known examples of bound 
entangled states show very low entanglement, therefore 
they will be more difficult to be checked experimentally. 

Now we present our first procedure, which checks en- 
tanglement in NPT bipartite states. The method we pro- 
pose can be thought of as a way of checking the Peres- 
Horodecki criterion. Assuming the knowledge of n (n < 



d — 1) of the parameters r k 



(fc = l,2,...,d'-l) 



in Eq. 1, we check the existence of a PPT state com- 
patible with the available information. The nonexistence 
of such a state witnesses the entanglement of the state 
of interest. This can be done by means of the following 
very simple SDP: 



determine g 



subject to < 



g>0 
Tr(g) = 1 
Q T >0 
Tr(gP k ) = 



(5) 



r fe , fc = l,2, 



This SDP is a feasibility program. g r stands for the par- 
tial transpose of g. When this program is infeasible, i.e., 
when there is no PPT state compatible with the avail- 
able data, we are certain that the unknown p is NPT 
and, therefore, entangled. 

It could happen that the state of interest has passed 
through some known decohcrence channel, which re- 
stricts the state's marginals to some known form. One 
example is the Werner states [15], which correspond to 
depolarized states whose marginals are maximally mixed. 
The program in Eq. 5 can be easily modified to include 
this additional information, which corresponds to further 
constraints in the SDP, namely: 

determine g 



subject to < 



t'i 



Tr(paf 



'') = (Tij), 



(4) 



g > 

Tr(g) = 

g r >0 

g a = Pa 

gb = Pb 

Tr{gP k ) 



1 



(6) 



r k , fc = l,2,...,n. 
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Programs in Eq. 5 and in Eq. 6 determine the projec- 
tion of the state of interest in a certain hyperplane (in 
the Hilbert-Schimidt space), and check the existence of 
the family of PPT states with the same projection. If 
there is no such state, it means that the measured state 
is NPT, and therefore entangled. This suffices to check 
entanglement in spaces 2x2 (qubit-qubit) and 2x3 
(qubit-qutrit) [2]. In larger spaces, this approach still 
works for the NPT states, but it will not detect entan- 
gled PPT states. Now we introduce a procedure that, 
in principle, can detect both free and bound entangled 
states. 

When the state of interest is in a space which allows for 
bound entanglement [7], we need an entanglement wit- 
ness to check if it is separable or not. If wc eliminate the 
constraint of positivity of the partial transpose in the pro- 
grams of Eqs. 5 and 6, namely , g r > 0, those programs 
return a state g, which can be PPT or not, compatible 
with the available data. We then build an optimal en- 
tanglement witness (W g ) to g, and use it to estimate the 
entanglement of p, i.e., Tr(W e p) ~ Tr{W g g). Remember 
that an EW is a Hermitian operator with non-negative 
expectation values on separable states, but which can 
have a negative expectation value on an entangled state, 
in this case, we say that the EW detects the entangled 
state. The optimal EW of a state yields the most nega- 
tive expectation value, when compared to any other EW 
of the same kind, therefore Tr(W e p) > Tr(W p p). Note 
that EWs can be chosen to correspond to different en- 
tanglement quantifiers [16]. The EWs in this work have 
the constraint Tr(W) — 1, and correspond to the ran- 
dom robustness [17, 18], which measures how resilient to 
white noise is the entanglement. 

We need an error bar to our entanglement estimate 
given by Tr(W e g). In order to do that, we rewrite Eq. 1 
as: 

n d 2 

Q = J2 r k p k+ E r i P r ( 7 ) 

k— 1 j—n+1 

The first summation corresponds to the known data. 
Let's call it Qknown- Of course, pknown = Qknown- The 
second summation is yielded either by the program in 
Eq. 5) or Eq. 6, and we call it Qunknown- Now we can 
write our entanglement estimate as: 

Tr(W eP ) = Tr{W e g) ± \Tr(W e g unknown )\. (8) 

The techniques we use to build the optimal EW are 
based on SDPs, and are described in [3]. 

In Fig. 1, we show how the method performs for 
two qubits, one qubit and one qutrit, and two qutrits. 
Each graph is built out of 10 4 random NPT entan- 
gled states. We plot the efficacy of entanglement de- 
tection (i.e., number of states detected as entangled di- 
vided by 10 4 ) in the sample of states, against the num- 
ber of measured non-local parameters (viz. Eq. 4). 
Every time the program in Eq. 5 or in Eq. 6 is in- 
feasible for a given state, it means that the measured 



data were sufficient to detect entanglement. In the case 
of two qutrits, we also test the EW approach of Eq. 
8. A state is considered successfully detected as en- 
tangled, when both Tr(W e g) + \Tr(W e g unk nown)\ and 
Tr(W e g) - \Tr(W e g unkn0 wn)\ are negative. About 70% 
of the states are detected as entangled, with an effort 
which is roughly half of a full state tomography. On one 
hand, the less entangled is the state, more information 
we need to infer its entanglement. Therefore, the graphs 
show 100% efficacy only when all the tomographic pa- 
rameters are measured. On the other hand, highly en- 
tangled states are detected with the knowledge of only a 
few non-local parameters. 

III. CHOOSING WHAT TO MEASURE IN THE 
LABORATORY 

In the last section, we have described the general idea 
behind detecting entanglement of unknown states, based 
on partial information. Our main goal was to show that 
the proposed data post-processing is effective. Now we 
want to discuss how our technique could be actually im- 
plemented with projective measurements. A sensible way 
to do this is by grouping the observables in the smallest 
number of maximal commuting classes. Commuting ob- 
servables share a common set of eigenvectors and, con- 
sequently, can be simultaneously measured (i.e., can be 
simultaneously diagonalized) . Therefore, such a scheme 
would yield the smallest number of complete projective 
measurements to be done. 

Let's fix the basis of observables. For a Hilbert space of 
dimension d s , we introduce the shift and clock operators, 
namely: 

d e -l 

X=J2\j + l){j\ (9) 

and 

Z = £exp(— %>01, (10) 
i=o 1 

where {\j}, j = 0, . . . , d s — 1} is an orthonormal basis. 
For dimension 2, these operators are the usual Pauli ma- 
trices. We also define: 

Y = XZ (11) 

and, for d > 2, 

V = XZ 2 . (12) 

Two-particle observables are now defined as tensor prod- 
ucts of powers of these operators. For two qubits, one 
qubit and one qutrit, and two qutrits, Table I shows the 
complete bases of observables, with a convenient labelling 
[19]. 
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FIG. 1: Fraction of success of entanglement detection against 
the number of measured non-local parameters (viz. Eq. 4) , 
for a sample of 10 4 random NPT states, using the approaches 
described in Eqs. 5 and 6, for two qubits, one qubit and one 
qutrit, and two qutrits. For two qutrits, we also show the 
results using the EW (Eq. 8). 
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TABLE I: Complete bases of observables in Hilbert spaces of 
dimensions 2(8)2, 2® 3 and 3<8>3. Each two-particle observable 
is the tensor product between one operator of the first line by 
one operator of the first column. 



Now we can group the observables of Table I in max- 
imally commuting classes. For two qubits, we have the 
5 classes [20] shown in Table II; for one qubit and one 
qutrit, there are 12 classes, as shown in Table III; and 
finally, for two qutrits, there are the 10 classes [21, 22] 
shown in Table IV. Note that, for two qubits and two 
qutrits, the classes are disjoint sets. In each case, the 
simultaneous eigenvectors of each class form a set of Mu- 
tually Unbiased Bases (MUB) [23, 24], in the sense that 
any two vectors of different bases have the same overlap's 
absolute value. In the case of the 2® 3 system, the classes 
are neither disjoint and nor minimal. With 5 observables 
in each class, the minimal number of classes, for a total 
of 35 distinct operators, should be 7. It is the case for 
the systems 2 ® 2 and 3 (g> 3, with 15 observables divided 
in 5 classes of 3 operators, and 80 observables divided in 
10 classes of 8 operators, respectively. It is conjectured 
that there is no informationally complete set (in the to- 
mographic sense) of MUBs for the 2 ® 3 system [25-27], 
and it is known that generalized Pauli matrices (which 
is our choice of observables) are not extensible to MUBs 
[28]. 

In Fig. 2, we repeat the calculations of section II, 
namely Eqs. 1, 5 and 8, but now using the MUB projec- 
tors for the two qubits and two qutrits. 

Though we do not know MUBs for the 2 eg) 3 system, 
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Ci = { 1 4 13 } 
C 2 = { 2 8 14 } 
C 3 = { 3 12 15 } 
C 4 = { 5 9 10 } 
C 5 = { 6 7 11 } 

TABLE II: Five maximally commuting classes of observables 
(viz. Table I) for two qubits. The common eigenvectors of 
each class form a set of MUBs. 
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TABLE III: Twelve maximally commuting classes of observ- 
ables (viz. Table I) for qubit®qutrit. 



we still want to do the minimal number of projective 
measurements in the laboratory. To make a complete to- 
mography, we need a set of 35 informationally complete 
projectors. Measuring in the basis of common eigenvec- 
tors of each of the 12 classes (Table III), the numbers 
of independent projective measurements extracted from 
each class are, respectively, 5, 4, 4, 4, 3, 2, 2, 2, 3, 2, 2, 
2. These 35 projectors, which are linearly independent 
in the Hilbert-Schimidt space, can be sorted in 7 sets of 
5, and re-orthonormalized in order to correspond to 7 
complete projective measurements (7 observables). The 
results of measurements in these two different bases are 
shown in Fig. 3. 
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TABLE IV: Ten maximally commuting classes of observables 
(viz. Table I) for two qutrits. The common eigenvectors of 
each class form a set of MUBs. 
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FIG. 2: Fraction of success of entanglement detection against 
the number of measured MUB projectors (viz. Eq. 1) , for 
a sample of 10 4 random NPT states, using the approaches 
described in Eq. 5 (Peres-Horodecki criterion), for two qubits 
(top), and Eqs. 5 and 8 (EW) for two qutrits (bottom). 
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FIG. 3: Fraction of success of entanglement detection against 
the number of optimal projectors (viz. Eq. 1 and Table III) , 
for a sample of 10 4 random NPT states, using the approach 
described in Eq. 5 (Peres-Horodecki criterion), for the 2 (g) 3 
system. 
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IV. ESTIMATED EW AND LOW-ENTANGLED 
STATES 

In sections II and III, we have discussed how we could 
detect the entanglement of unknown NPT states, based 
on partial information. In particular, the method we have 
proposed to check the Peres-Horodecki criterion (Eq. 5) 
is rigorous, yields an exact answer and, according to our 
numerical tests, performs nicely for the systems 2 <g> 2, 
2 (g) 3 and 3 3, as shown in Figs. 2 and 3. On the 
other hand, our proposed EW estimate (Eq. 8) needs 
to be better understood. Figs. 1 and 2 are numerical 
evidence for the correctness of Eq. 8, but we are lacking 
a rigorous proof for the exact expression of our error 
bar. In this section, we will study the performance of 
Eq. 8 in 3 particular states, being one highly entangled 
two-qutrit Werner state [15] and two very low entangled 
states, being one of them also a two-qutrit Werner state, 
and the other one a two-qutrit bound entangled state 
[29]. It will add further evidence of the correctness of 
Eq. 8, and will show that our proposed error bar is too 
big, i.e., it seems that Tr{W e g) is a very good upper 
bound to Tr(W e p), much better than we expected, and 
there must be a tighter error bar, but we couldn't devise 
it yet. Note that Tr(W e p) is certainly an upper bound 
to Tr{W p p). 

The two-qudit (for our purposes d = 3) Werner states 
[15] can be written as: 



Pw 



i d + m 

d 2 + d/3 ' 



with — 1 < (3 < 1. pw is separable for (3 > 
swap operator for two qudits, 



(13) 



F d is a 



(14) 



The two-qutrit bound entangled state we use is picked 
up from the following family of states [29]: 
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PH = ^i)(^\ + ^a + + ^a_, 2<A<5, (15) 
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where 



X i,3=0 



(16) 



is the density matrix for the maximally entangled state, 
and 

a + = i(|01)(01| + |12)(12| + |20)<20|), (17) 



«7_ = -(|10)(10| + |21>(21| + |02)<02|) 



(18) 



FIG. 4: Estimated EW with its error bar (Eq. 8) for 3 par- 
ticular two-qutrit states: (top) Werner state with /3 = — 1, 
(middle) Werner state with (5 = —0.37, (bottom) Horodecki 
bound entangled state with A = 3.9. In each case, the exact 
value for Tr(W p p) corresponds to the mark 80 in the abscissa. 



are two separable states. With these definitions, the char- 
acter of ph changes with A according to 

separable, 2 < A < 3, 

Ph = < bound entangled, 3 < A < 4, (19) 
free entangled, 4 < A < 5. 
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In Fig. 4, we see the results yielded by Eq. 8 applied 
to Werner states (Eq. 13) with = -1 (Tr(W pw p w ) = 
-1/3) and (3 = -0.37 (Tr(W Pw p w ) = -0.014), and 
to the Horodecki bound entangled state (Eq. 15) with 
A = 3.9 (Tr(W PH p H ) = -0.024). As we have mentioned 
before, a state is considered detected as entangled when 
the error bar resides entirely in the entangled region. We 
see that the highly entangled Werner state is detected 
with just 11 measurements. On the other hand, the two 
low entangled states are detected after the 70*^ measure- 
ment. The Peres-Horodecki criterion (Eq. 5) detects the 
low entangled Werner state in the 32 nd measurement but, 
of course, it is not applicable to the PPT state. Ignoring 
the error bar, note that the estimated EW never super- 
estimated the entanglement; it yielded the exact results 
for the Werner states after the 17 th measurement, for the 
most entangled state, and after the 32 nd measurement, 
for the low entangled state; and, finally, it detected the 
bound entangled state after 10 measurements. 

V. CONCLUSION 

We discussed data post-processing strategies to char- 
acterize entanglement of unknown mixed states, based on 
partial knowledge of the state. The method is guaran- 
teed to work, for it converges to a full state tomography. 
We applied our method in systems of dimension 2 ® 2, 
2(g) 3 and 3(g) 3. Our numerical investigations showed that 
entanglement can be detected with a cost which is much 
lower than full state tomography, when the entanglement 



is not very small. For low entangled states, including 
PPT ones, we presented a method to construct entangle- 
ment witnesses (EW). The EWs have an error bar that 
monotonically diminishes with the increase of informa- 
tion about the state. Our tests suggest that the error 
bar is too big, for ignoring it, the entanglement estimate 
yielded by the EW is always a lower bound to the true 
entanglement. Therefore we believe that a tighter error 
bar could be calculated, but we weren't able to prove it 
yet. 

We also discussed the choice of observables to be mea- 
sured in the laboratory. In particular, we noted that the 
choice is not obvious in the case of the 2(g3 system, when 
one is willing to measure the smallest set of informa- 
tionally complete projectors. Nevertheless, we offered a 
method to construct these minimal informationally com- 
plete sets, in the case of projective measurements. 

The application of our approach to multipartite sys- 
tems is straightforward, at the level of the formalism. As 
a matter of fact, we performed some tests on NPT states 
of three qubits, obtaining results similar to the ones we 
presented for the bipartite systems. 

It would be interesting to investigate how our approach 
could be adapted to study other properties, as the purity 
of a state, for example. In this case, quadratic or higher 
degree non-linear programs would be necessary. 
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